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After fixing the Maximal Abelian gauge in SU (2) lattice gauge theory we decompose the nonabelian gauge field 
into the so called monopole field and the modified nonabelian field with monopoles removed. We then calcnlate 
respective static potentials and find that the potential due to the modified nonabelian field is nonconfining while, 
as is well known, the monopole field potential is linear. Furthermore, we show that the sum of these potentials 
approximates the nonabelian static potential with 5% or higher precision at all distances considered. We conclude 
that at large distances the monopole field potential describes the classical energy of the hadronic string while 
the modified nonabelian field potential describes the string fluctuations. Similar decomposition was observed to 
work for the adjoint static potential. A check was also made of the center projection in the direct center gauge. 
Two static potentials, determined by projected Z 2 and by modified nonabelian field without Z 2 component were 
calculated. It was found that their sum is a substantially worse approximation of the SU{2) static potential than 
that found in the monopole case. It is further demonstrated that similar decomposition can be made for the flux 
tube action/energy density. 


1. Introduction 

It is well known in™ that after Abelian 
projection in the Maximal Abelian gauge (MAG) 
|5I6| the abelian string tension, calculated from 
the Abelian static potential, is very close to 
the nonabelian string tension and corresponding 
Coulomb term coefficient is about 1/3 of that in 
the nonabelian static potential. The former ob¬ 
servation as many others supports the Abelian 
dominance (for review see e.g. |7]). It was fur¬ 
ther discovered that the monopole static 

potential also has a string tension close to the 
nonabelian one and small Coulomb term coeffi¬ 
cient. These observations are in agreement with 
conjecture that monopole degrees of freedom are 
responsible for confinement cni. 

The role of the other, nonrelated to monopoles, 
degrees of freedom was investigated to a lesser ex- 
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tent so far. We believe that studying the prop¬ 
erties of the observables constructed out of these 
nonmonopole degrees of freedom will help to un¬ 
derstand the properties of the infrared effective 
action. 

It is then interesting to see what kind of static 
potential one gets if one switches off the monopole 
contribution to the gauge field, i.e. if only off- 
diagonal gluons and the so called photon part of 
the abelian gluon field are left interacting with 
static quarks. 

Previously computations of this kind were 
made in inna, where it was shown that the 
topological charge, chiral condensate and effects 
of chiral symmetry breaking in quenched light 
hadron spectrum disappear after removal of the 
monopole contribution from the relevant opera¬ 
tors. Similar computations were made within the 
scope of the Z 2 projection studies m It was 
shown that the modified gauge field with pro¬ 
jected center vortices (P-vortices) removed pro- 
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duces Wilson loops without area law, i.e. lack¬ 
ing confinement property. In fact we want to do 
a similar thing with monopoles. But we go one 
step further - we consider the sum of the static 
potentials Vmod(R) and VmoniR), obtained from 
the Wilson loops of the modified gauge field (with 
monopole contribution removed) and from the 
Wilson loops of the monopole gauge field alone, 
respectively. We discover that Vmod{R) can be 
well fitted by pure Coulomb term and the sum 
VmodiR) + Vmon{R) is a good approximation of 
the nonabelian static potential V{R) at all dis¬ 
tances considered in our measurements. Further¬ 
more, we find that a similar approximate decom¬ 
position holds for the adjoint static potential. We 
also check the decomposition of the static po¬ 
tential induced by the center projection. The 
approximate decomposition for the static poten¬ 
tial suggests the possibility of such decomposition 
for the flux tube action/energy density. We also 
check this possibility. 

The paper is organized as follows. In the next 
section we introduce the necessary definitions and 
describe details of our computations. In section 
3 results for the various static potentials are pre¬ 
sented. Section 4 is devoted to the action/energy 
density results. Finally, we conclude in section 5. 

2. Definitions and simulation details 

We study SU(2) lattice gauge theory with Wil¬ 
son action. The abelian projection is made af¬ 
ter fixing to the MAG. The abelian projection 
means coset decomposition of the nonabelian lat¬ 
tice gauge field t/(s,/j,) into abelian field u{s,^) 
and coset field C{s,fi): 

U{s,fj,) =C{s,n)u{s,fi). ( 1 ) 

The Abelian gauge field can be decomposed into a 
monopole (singular) part Umon{s,^) and photon 
(regular) part Uphis,fj,): 

^(^5 M) — (2) 

or for the corresponding angles: 

0{s, H) = 0monis, ^J.) + 0phis, , (3) 

where 0{s, ^),0mon,phis, G (—7r,7r] are defined 

by relations 

u{s,fj.) = (4) 


w^on.p/.(s,M) = (5) 

0mon{s,^) satisfies the equation |14j : 

^L'^i/0mon{,^: f^') ^fi^i/0mon(^^t — ( 6 ) 

2'Kd^m{s, vn) , 

where dp (d^) are lattice forward (backward) 
derivatives. The Dirac plaquette variable 
m(s, fdv) G Z is determined by decomposition of 
the abelian plaquette angle 0{s, iiu) = dp0{s, v) — 
d^0{s, into regular and singular parts: 

0{s, ^ly) = 0{s, fiv) + 2Trm(s, fiv), (7) 

0{s,fj,iy) G (—7r,7r] 

Equation m has solution 

0Tnonis, n) = -27r ^ D{s - s')d^m{s', vy) , (8) 

s' 

where D{s) is the lattice inverse Laplacian. This 
solution satisfies the Landau gauge condition 
dp0mon{s, fi) = 0. The monopole gauge field 
0mon{s, /r) defined in eq.® reproduces more than 
95% of monopoles on our lattices and it repro¬ 
duces all monopoles in the infinite volume limit. 
This explains its name. We calculate the usual 
Wilson loops: 

W{C) = ^Tt(i[u{1)Y (9) 

Viec / 

as well as monopole Wilson loops: 

W^oniC) = ^Tr Umonil)^ , (10) 

and the nonabelian Wilson loops with monopole 
contribution removed: 

kE„^od(C) = pr|^n^(0^ ’ 

where the modified nonabelian gauge field is de¬ 
fined as 

U{s,fj,) = C{s,fi)uph{s,fJ,) . (12) 

Note that Uph{s,^) is abelian projection of 
U{s, fj,) and it has very few monopoles in the finite 
volume and no monopoles in the infinite volume. 
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Fixing MAG leaves unbroken U{1) gauge sym¬ 
metry. The general form of this transformation 
is 

9'{s, /r) = 9{s, /r) -I- df^uj{s) + 27rn(s, fi), (13) 

where 9'{s, ^),uj(s) S (—tt, tt], n(s,fi) = 
0, ±1. Respectively, the Dirac plaquette variable 
changes as 

to '( s , V, = m(s, V, ^)-\-dvn[s^ ^)—d^n{s, v) .(14) 

Thus the monopole field eq. 0 depends on the 
choice of the U{1) gauge: 

^mon(s.M) = 9mon{s, fl) + 2Trn{s, n) - 6{s, fl) , (15) 
where 

i5(s,/i) = — 2'k'^^D{s — s')dfj,d^n{s',v) (16) 

s' 

- TT—XI 

sites g 

The monopole Wilson loop Wmon{C) is invariant 
under the change of the monopole field 113 due 
to the fact that ~ 

for WmodiC) because of its nonabelian character. 
Indeed, the transformation for U(s,^) is 

U'{s,fi) = n{s)U{s,fi)n{s + fj.)^A{s,n), (18) 

where 

n(s) = diag{e(-“("», , (19) 

A(s, fi) = diag{e-*'^("’'^), } . (20) 

This transformation is not a gauge transforma¬ 
tion and thus the Wilson loop Wmod depends on 
the choice of the gauge used in the definition of 
the monopole field eq. ©■ This problem is solved 
by using in m the Dirac plaquettes m(s, ^v) in 
some particular gauge. We choose the Landau 
gauge defined by the gauge condition: 

maxY^ cos{9'{s, y)). (21) 

U) 

S,^ 

Another possible gauge condition would be the 
minimization of the number of the Dirac plaque¬ 
ttes: 

min (m'(s,/r, i^))^ . (22) 


Up to Gribov copies both conditions fix config¬ 
uration of Dirac plaquettes m{s,^v) completely 
and thus fix 9mon{s, fi). In general, results for 
Wmod can be different for different gauge con¬ 
ditions which fix the Dirac plaquette configura¬ 
tion. But for two gauge conditions introduced 
above we may hope that results are similar be¬ 
cause fixing t/(l) Landau gauge also strongly re¬ 
duces the number of Dirac plaquettes and thus we 
expect that for given lattice configuration num¬ 
ber of Dirac plaquettes in these two gauges are 
close to each other. 

We calculated RxT rectangular Wilson loops 
W{R,T), Wmon{R,T) and Wmod{R,T). To ex¬ 
tract the nonabelian static potential V{R) link 
integration m and smearing m have been em¬ 
ployed. Smearing was used in computations of 
the modified field potential Vmod(,R)- Gomputa- 
tions were done at /3 = 2.5 on 24"^ lattices us¬ 
ing 100 statistically independent configurations. 
To fix MAG 10 randomly generated gauge copies 
fixed by simulating annealing algorithm were 
used. 

3. Static potential decomposition 

In Fig.nileft) we show the monopole Vmon{R) 
and the modified field Vmod{R) potentials. We 
find that Vmon (R) is linear at large distances and 
has small curvature at small distances, as was ob¬ 
served many times before. Our result for Vmod{R) 
is the first result for this potential. It can be 
seen from Fig. n(left) that this potential is of the 
Goulombic type. Indeed it can be very well fitted 
by Y^^-amodlR with Umod = 0.274(9). The fit¬ 
ting curve is shown in Fig. n^left). Thus, remov¬ 
ing the monopole contribution from the Wilson 
loop operator leaves us with Wilson loop which 
has no area law behavior, i.e. confinement prop¬ 
erty is lost. This result is similar to that obtained 
in m after removing P-vortices. 

In Fig. n^left) we also compare the sum 
Vmon(R) + VmodiR) with the nonabelian poten¬ 
tial V (R). One can see that the nonabelian static 
potential is well approximated by this sum, i.e. 

ViR) « VmoniR) + VmodiR) ■ ( 23 ) 

This is our main result. On the quantitative level. 
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Figure 1. (left) Comparison of the nonabelian potential V{R) (squares) and the sum Vmod{R) + Vmon{R) 
(triangles). Vmod{R) (diamonds) and Vmon{R) (circles) are also depicted, (right) Vmod after U(l) gauge 
fixing (diamonds) and without U(l) gauge fixing (triangles). 


we find that 


5V{R) 

V{R) 


< 0.05, 


The approximate relation implies that for 
large T : 

{W{R,T)) « K {W^on{R,T)) (WmodiR,T)) (25) 


where 6V{R) = V{R) - {V,^on{R) + V^odiR))■ 
This observation can be formulated in the follow¬ 
ing way: the potential for the static sources, inter¬ 
acting with the nonabelian gauge field can 

be approximated by the sum of the potential for 
the sources, interacting only with the monopole 
field Umonis, /r), and potential for the sources, in¬ 
teracting only with the modified field U{s, ^). All 
calculated potentials were fitted by the usual lin¬ 
ear plus Coulomb functions. Results for the fit 
parameters are presented in Tabled 

In Fig. fright) we compare Vmod{R) calculated 
with and without additional U(l) gauge fixing de¬ 
fined by eq. m- One can see that without U(l) 
gauge fixing Vmod(,R) is substantially lower than 
it is after U(l) gauge fixing. Fitting with con¬ 
stant plus Coulomb fitting function we found the 
Coulomb coefficient two times smaller. Thus the 
approximate decomposition of the potential ra 
does not hold without this gauge fixing. 


Indeed we found that for A, T > 3 our data for 
unsmeared {W{R,T)) can be fitted by the right 
hand side of eq. (1^ (with unsmeared Wmon and 
Wmod) with K = 0.886(9). 

Next we want to address the adjoint potential 
decomposition. We calculated nonabelian adjoint 
Wilson loop WadjiR^T), charge two monopole 
Wilson loop 

W^ou,2i.R.T) = Wl,^{R,T) (26) 

and adjoint Wilson loop Wmod,adjiR,T) for the 
modified nonabelian field U{s,^). From these 
Wilson loops the respective potentials Vadj{R), 
VmonMR) and Vmod,adjiR) extracted. In 
Fig. |2Ileft) we compare the sum Vmon, 2 {R) + 
Vmod,adj{R) ■^ith the adjoint potential Vadj{R)- 
As in the fundamental case we see approximate 
decomposition: 

^adj^R') ~ ^mon,2{R') “f ^mod^adj ^R') ■ 


(27) 












Table 1 

Parameters of the potentials obtained by fits with function Vq — a/i? + aR. 
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aa^ 

a 

aVo 

aV{R) 

0.0339(4) 

0.286(4) 

0.560(2) 


0.0310(3) 

0.264(6) 

0.577(4) 

<3'Tdnon(-f^) 

0.0311(3) 

-0.018(4) 

-0.006(3) 

^^mod (^) 

-0.0002(2) 

0.280(10) 

0.583(5) 

aVz^ (R) + aVmod,Z2{R) 

0.0224(15) 

0.334(20) 

0.668(15) 

aVz^iR) 

0.0249(6) 

-0.024(7) 

-0.003(4) 

^^7nod,Z2 (^) 

-0.029(6) 

0.383(10) 

0.686(6) 


The potential Vmod,adj[R) looks purely Coulom- 
bic. A Coulomb term coefficient amod,adj = 
0.69(4) was found by fitting the data at i? > 2. 
Vmon, 2 (R) is linear with small curvature as was 
observed before in |S]. Our result supports 
the conjecture IT7I^ that the abelian charge two 
potential should be considered as the abelian pro¬ 
jection for the adjoint potential. 

As we noted above, in Ref. m it was shown 
that in the central gauge after removal of P- 
vortices the Wilson loop looses confinement prop¬ 
erty. In the MAG we found that removal of 
monopoles also leaves the Wilson loop without 
confinement property. Moreover, we found that 
the approximate decomposition holds. It is 
then interesting to check similar decomposition 
for the central gauge. We made computations in 
the Direct Central (DC) gauge at /3 = 2.5 us¬ 
ing half of the set of configurations which were 
used for computations in the MAG. Results are 
presented in Fig. I^^right). From comparison of 
Fig. ^left) and Fig. [fright) and from compari¬ 
son of the respective fitting parameters, presented 
in Tabled one can see that approximate decom¬ 
position works substantially better in MAG. This 
is not unexpected for large distances, since it is 
known that in the center gauges the P-vortex 
string tension is smaller than the monopole string 
tension in MAG. 

Note that decomposition similar to eq. m is 
impossible for central projection since the charge 
two central projected Wilson loop is identity. 


4. Flux tube profile 

In this section we present our results for the 
flux tube action/energy density decomposition. 

Taking derivative with respect to (3 on the left 
and right hand sides of eq. (1^ and dividing by 
respective expressions we obtain (for T —> oo) 

{SW{R,T)} _ _ f {SWmon{R,T)) _ \ 

W{R,T) W^oniR,T) 


({SWrnodiR^ _ , , 

V Wmod{R,T) ^ ’ 


(28) 


where S is the nonabelian action averaged over 
the lattice. It is then interesting to check whether 
such approximate decomposition is valid also lo¬ 
cally, i.e. whether the following approximate re¬ 
lation holds: 




i^Lnis, t^))w + {F^odis, ’^))w , (29) 


where 


{Ff{s,fj.,i2))w = 


{WfFHs,fd,i2)) 

(Wf) 




the lattice field strength squared is defined as 

^F^{s,^i,ly) = 1-^Tr U{s,^i,ly), (30) 

and F^„^(s,/r, j^) are defined anal- 
ogously with respective plaquettes instead of 


^In this derivation we take into account that the Faddeev- 
Popov determinant has no explicit dependence on 
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Figure 2. (left) The adjoint nonabelian potential Vadj{R) (squares) is compared with the sum 
Vmod,adj{R) + Vmon, 2 {R) (triangles). VmodMjW (diamonds) and Vmon, 2 {R) (circles) are also shown, 
(right) Comparison of the nonabelian potential V{R) (squares) and the sum Vmod,Z 2 {R) + Vz^iR) 
(triangles). Also depicted are Vmod,Z 2 {R) (diamonds) and Vz 2 {R) (circles). 


U{s,fj,,v). In Fig. 13 and Fig. 0 we show various 
longitudinal action and energy densities defined 
as follows: 

AL{r±,z) = - ((F^(s, 0, 3))u/ + (F^(s, 1, 2))w) , 

EL{r±,z) = - ({F^{s,0,3))w - {F^is,l,2))w) , 

where r_L is the distance from the quark-anti¬ 
quark axes, 2 is coordinate along this axes. The 
distance between sources is R/a = 8. One can 
see from these figures that relation indeed 
holds with rather good precision. Similar results 
were obtained for the transverse action and en¬ 
ergy densities. In the computations of the ac¬ 
tion/energy density Wilson loops with T/a = 4 
were used for the nonabelian and modified cases 
and with Tja = 6 for the monopole case. 

5. Conclusions and discussion 

Using abelian projection after fixing MAG we 
shown that the static potential V(R) can be ap¬ 
proximately decomposed into two components, 


see eq. (P|l. Our preliminary results at smaller 
lattice spacing (/3 = 2.6) indicate that in the 
continuum limit the agreement between left and 
right hand sides in (PI improves. One term in 
(P|. VmoniR), is due to the monopole gauge field 
(0 contribution, another one, Vmod{R), is due to 
the contribution of the modified nonabelian gauge 
field U{s,fj,) with monopoles removed. Compar¬ 
ison of the parameters of these potentials, given 
in Table ^ with the effective string model pre¬ 
dictions suggests the following interpretation: at 
large distances Vmon{R) is the classical energy of 
the flux tube and Vmod{R) is the flux tube fluctu¬ 
ations energy. This conclusion is supported by 
our results for the adjoint static potential, see 
Fig.01eft). 

It was checked whether a similar decomposition 
of V(R) holds for the center projection after fixing 
to DC gauge. We calculated projected Z 2 poten¬ 
tial Vz 2 (R) and modified nonabelian field poten¬ 
tial Vmod,Z 2 (R)- The latter was previously intro¬ 
duced and measured in m We found that de¬ 
composition in this case holds with substantially 
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z r_L 

Figure 3. The longitudinal action densities: nonabelian (squares), monopole (diamonds), and sum of the 
monopole and modified (circles), (left): as functions of the coordinate z along the QQ axes; (right): as 
functions of the distance r± from this axis. 
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Figure 4. Same as in Fig. 0 but for the longitudinal energy densities. 






























less precision than in the monopole case. This 
can be seen from comparison of Figs. mieft) and 
|2Iright), and from Tabled Both the string ten¬ 
sion and the Coulomb coefficient obtained from 
the hts of the sum Vz^iR) + Vmod,Z 2 {R) are 
in worse agreement with respective parameters 
of the nonabelian potential than those obtained 
from the fits of the sum Vmon{R) + Vmod{R)- 

The decompositions eq. and eq. are 
similar to decomposition in the compact f7(l) 
model m into the photon and monopole com¬ 
ponents of the static potential. In that model the 
action can be respectively decomposed into the 
monopole and photon parts without interaction 
term. Contrary to the compact U{1) theory such 
term is unavoidable in the infrared effective ac¬ 
tion of the SU{2) gluodynamics since the charge 
two component of the off-diagonal held C{s, /i) 
provides the string breaking of the charge two 
monopole potential Vmon ,2 Hi The relations 
eq. (1^ and eq. (E7I) are hinting that this inter¬ 
action term should be weak. 

Our results for the action and energy densities 
in the quark-antiquark system, presented in Fig.|31 
and Fig. 21 show that the approximate decompo¬ 
sition holds also for these quantities. 
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